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Abstract 


New analyses of earlier ALMA observations of oxygen-rich AGB star EP Aquarii are presented, which complete a 
previously published analysis and offer a different interpretation of the morpho-kinematics of the circumstellar 
envelope. The birth of the equatorial density enhancement (EDE) is shown to occur very close to the star where 
evidence for rotation has been obtained. Close to the star and where outflows have been observed: their interaction 
with the gas of the nascent EDE is seen to play an important role in the development of the wind and the evolution 
of its radial velocity from 8 to 10kms~' on the polar symmetry axis to —2 kms '' at the equator. It implies 
complex morpho-kinematics: making reliable interpretations with reasonable confidence is difficult. In particular, it 
questions an earlier interpretation implying the presence of a white dwarf companion orbiting the star at an angular 
distance of ~0”4 from its center. It proposes instead an interpretation in terms of a standard mass ejection 
associated with a shock wave leaving a void of emission in its wake. High Doppler velocity wings are seen to 
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consist of two components, the upper velocity end of the global wind, reaching above +12 km s ^ !, and an effective 
line broadening, confined within 200 mas from the center of the star, reaching above 4-20 km s! and interpreted as 
caused by the pattern of shock waves resulting from the interaction between stellar pulsation and convective cell 


partition. 
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1. Introduction 


When instruments providing observations of high sensitivity 
and angular and spectral resolutions became available, mainly 
the Very Large Telescope (VLT) and the Atacama Large 
Millimeter/submillimeter Array (ALMA), studies of the circum- 
stellar envelopes (CSEs) of Asymptotic Giant Branch (AGB) 
stars progressed rapidly. Yet, many questions are still unan- 
swered concerning the mechanism that governs the generation 
and acceleration of the wind, in particular for oxygen-rich low 
mass-loss rate stars. In contrast with carbon-rich stars, these have 
a dusty envelope that is more transparent in the visible and near- 
infrared (near-IR), giving a better view of the innermost dust- 
forming atmospheric layers (Hófner & Freytag 2019). Their 
study supports a scenario where the radiation pressure, which 
accelerates the outflows, is caused by photon scattering on 
highly transparent iron-free grains, such as corundum (ALO), 
which form a thin gravitationally bound dust layer close to the 
stellar photosphere; beyond it, silicate dust condenses on top of 
corundum cores, speeding up grain growth to the critical size. 

Recent high-resolution imaging of nearby AGB stars at 
visible and infrared wavelengths has revealed complex, non- 
spherical distributions of gas and dust in the close circumstellar 
environment, with changes in morphology and grain sizes 
occurring over the course of weeks or months. At millimeter 


wavelengths complexity is the main feature that characterizes 
the morpho-kinematics of the CSE, with important anisotropy 
and inhomogeneity. Such distribution of atmospheric gas 
emerges naturally in 3D hydrodynamical models (Hófner & 
Freytag 2019) as a consequence of large-scale convective flows 
below the photosphere and the resulting network of atmo- 
spheric shock waves, the dynamical patterns in the gas being 
imprinted on the dust in the close stellar environment. Global 
AGB star models are characterized by giant convection cells, 
which can span over a steradian, and have a lifetime of many 
years. Together with radial pulsations, with typical periods at a 
few-month scale, they generate waves of various frequencies 
and spatial scales, which quickly develop into shock waves 
giving rise to ballistic gas motions peaking around 2 stellar 
radii. As the shock waves interact and merge, they produce 
large-scale regions of enhanced densities in their wakes. 

A consequence of the resulting complexity is the difficulty to 
reconstruct reliably the observed morpho-kinematics of the CSE in 
both velocity and space, driving its study into proceeding 
somewhat by trials and errors: successive observations and 
analyses bring new information that helps refining the picture. 
EP Aquarii (EP Aqr) provides such an example with four recent 
publications by Homan et al. (2018), Hoai et al. (2019), Tuan-Anh 
et al. (2019) and Homan et al. (20202). It is an M-type, semi- 


Research in Astronomy and Astrophysics, 24:035009 (23pp), 2024 March 


regular variable star (spectral type M8III), with a pulsation period 
usually quoted as 55days but more like 110-20 days 
(Eggen 1973; Tabur et al. 2009), located at 119 + 6 pc from the 
solar system (van Leeuwen 2007; Gaia Collaboration 2018); it 
shows no technetium in its spectrum (Lebzelter & Hron 1999) and 
has an effective temperature of — 3240 K, a mass of ~1.7 Mo, a 
radius of ~0.77 au (Dumm & Schild 1998) and a mass-loss rate of 
«(1.6 x: 0.4) 107 M yr ' (Hoai et al. 2019). 

The resulting picture that emerges from the four more recent 
studies of the CSE is two-component morpho-kinematics 
producing a Doppler velocity spectrum made of a narrow central 
peak at star systemic velocity and a broader plateau extending 
over ~ + 10-11 kms '. The central peak is associated with an 
oblate, and the broad plateau with a prolate, volume of gas, both 
expanding radially and sharing a same axis that projects ~20° 
west of north and is inclined by only —10? with respect to the 
line of sight. The oblate volume of gas, referred to as equatorial 
density enhancement (EDE), or occasionally torus or disk, is 
nearly face on, causing its radial expansion velocity to be poorly 
defined, —-2kms '. It hosts a complex pattern of arcs of 
enhanced CO(2-1) emission that is interpreted by Homan et al. 
(2018, 20202) in terms of a spiral produced by a white dwarf 
companion, with a mass of 0.65-0.80 Mo, orbiting EP Aqr at a 
distance of ~0”4 and currently located west of it. The prolate 
volume of gas, usually referred to as polar or bi-conical outflows, 
displays maximal velocity along the line of poles but clear 
density depressions near the poles, with the maximum being 
reached at stellar latitudes between ~ +45° and — +60° while 
the transition between oblate and prolate components occurs 
progressively near stellar latitudes of ~ +20°. 

Evidence for rotation within some 100—200 mas from the 
center of the star, with a rotation velocity at km s! scale, has 
been obtained from the observation of the emissions of the 
COQ-1), SO»(1656.10 = 17543) and SO5(455 zt 313) lines 
(Homan et al. 2018; Tuan-Anh et al. 2019; Homan et al. 
2020a). Observations of the dust induced polarization of the 
550-750 nm emission observed at the VLT using SPHERE- 
ZIMPOL (Homan et al. 20202) give evidence for a shell of dust 
surrounding the star and having an inner radius slightly smaller 
than 100 mas. A study of dust emission in the infrared using the 
Short Wavelength Spectrometer aboard Infrared Space Obser- 
vatory (ISO, Heras & Hony 2005) has given evidence for very 
low ratios of both gas-to-dust mass and aluminum oxide to 
silicates relative abundance. 

There exists no well-established interpretation of the observed 
morphology. EDE's are usually described as the result of the 
gravitational attraction of the stellar wind in the orbital plane of a 
companion and the polar outflows are often interpreted as 
episodically produced by the rapid rotation of an accretion disk 
surrounding this companion and reach very high velocities. Such 
are the cases of HD 101584 (Olofsson et al. 2019) with highly 
collimated outflows reaching velocities of ~ + 140 km s Lofa! 
Gruis (Doan et al. 2017; Homan et al. 2020b; Doan et al. 2020) 
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with polar outflows consisting of thin hourglass-shaped bubbles 
reaching a radial expansion velocity of ~60kms_', of V Hya 
(Sahai et al. 2022) with high velocity lumpy jets of velocity 
~175 kms '!, and of R Aqr (Liimets et al. 2018; Melnikov 
et al. 2018; Bujarrabal et al. 2021) with outflows made of curved 
central jets with velocity reaching ~240kms_' surrounded by 
hourglass bubbles in ballistic expansion with velocity of 
~55kms_'. Much closer to the morpho-kinematics of the EP 
Aqr CSE is that of AGB star RS Cnc, which however displays a 
clear technetium signal, but is too boreal to be observed by 
ALMA and has been studied in detail using NOEMA (Winters 
et al. 2022). It has two-component morpho-kinematics, quantita- 
tively similar to EP Agr, with an axis making an angle of ~30° 
with the line of sight and radial expansion velocities reaching 
3-4 kms ! in the EDE and 8-9 kms ' in the bipolar outflows. 

Neither EP Aqr nor RS Cnc has been shown to host a 
companion (but in both cases the presence of an undetected 
companion cannot be excluded) and no other known AGB star 
displays close enough morpho-kinematics to be claimed to 
serve as a model for their CSE. The question of their genesis 
remains therefore open and, in the case of EP Aqr, the picture 
proposed by Homan et al. (20202), suggesting the presence of 
an undetected white dwarf companion playing a major role in 
the formation and structure of the EDE, leaves several 
questions unanswered that call for further study: 

(i) High Doppler velocity wings, reaching ~ + 20km E 
were interpreted by Tuan-Anh et al. (2019) as two narrow polar 
structures, referred to as jets, launched from less than 25 au 
away from the star and building up between ~20 and ~100 au 
to maximal velocity. However, subsequent analyses of similar 
high Doppler velocity wings observed in several other stars, as, 
for example o Ceti (Hoai et al. 2020; Nhung et al. 2022), 
favored interpretations in terms of shocks produced close to the 
star by pulsations and convective cells, shedding strong doubts 
on the jet interpretation proposed for EP Aqr. It has therefore 
become imperative to reconsider the case. 

(ii) The formation of the bipolar outflows remains unex- 
plained. It may invite interpretations in terms of magnetic fields 
(García-Segura et al. 2005; Winters et al. 2022), or in terms of 
interacting winds in a binary (Nordhaus & Blackman 2006; 
Castellanos-Ramírez et al. 2021), and more specifically in 
terms of an additional companion orbiting EP Aqr at a distance 
of ~10au (Homan et al. 2020a), but none of these can be 
solidly justified. The latter authors present a number of relevant 
suggestive arguments in relation to the observation of a 
complex pattern of arcs in the emission of the polar outflows. 

(iii) The formation of the EDE remains equally unexplained. 
The interpretation proposed by Homan et al. (20202), in terms 
of focusing by an assumed white dwarf companion, is in 
principle compelling but raises unanswered questions and is 
likely to face inconsistencies when details of the observed 
morpho-kinematics are being considered. 
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Figure 1. (a) intensity map of CO(2-1) emission; (b) intensity map of SiO(5-4) emission; (c) Doppler velocity spectra of both CO(2-1) and SiO(5-4) line emissions; (d) 


schematic of the redshifts resulting from the tilted axi-symmetric geometry. The color scales are in units of Jy arcsec ? km s7". 
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Table 1 
Summary of Band 6 Observations and Data Reduction of the CO(2-1) (230.55 GHz) and SiO(5-4) (217.12 GHz) line Emissions 
Reference Epoch Line Beam Size Beam PA Channel Width Time Noise rms Maximal Recoverable 
(mas?) (°) (km s7’) on source (mJy beam!) Scale (”) 
Homan et al. (2018) 2016 COQ-1) 170 x 150 —34 00.33 70 minutes (ACA) 1.25 15 
SiO(5-4) 230 minutes (TM1,2) 

Combined data 2016, COQ-1) 87 x 79 +15 ~0.33 es 0.65 15 

Homan et al. (20202) 2019 SiO(5-4) 85 x 79 +7 

TMO data 2019 CO(2-1) 32 x 28 +13 0.95 99 minutes 1.0 0.5 

present work SiO(5-4) 32 x 30 +22 0.34 1.34 0.5 


(iv) A clear depression of SiO(5-4) line emission is observed 
in the western blueshifted octant and interpreted by Homan 
et al. (2020a) as the result of molecular dissociation produced 
by the ultraviolet (UV) radiation of the white dwarf companion. 
However, to the extent that the companion can be naively 
expected to sit in the median plane of the disk that it has 
generated, such dissociation should extend to the redshifted 
hemisphere as well and the observed morphology of the 
depression is difficult to understand in such terms. 

The aim of the present article is therefore to shed new light 
on these issues, by presenting new additional analyses of the 
emissions of the CO(2-1) and SiO(5-4) lines, in particular in 
the close neighborhood of the star. As a brief reminder of the 
main properties of the CSE of EP Aqr, we display in Figure 1 
intensity maps of the CO(2-1) and SiO(5-4) line emissions, 
together with Doppler velocity spectra and a schematic 
description of the axi-symmetric geometry. 


2. Observations and Data Reduction 


We use archival ALMA observations of the CO(2-1) and SiO 
(5-4) line emissions of EP Aqr that were performed in 2016 and 
2019 under projects ADS/JAO.ALMA#2016.1.00057.S and 
ADS/JAO.ALMAZ22018.1.00750.S, respectively (PI W. 
Homan). Observations and data reduction have been described 
in much detail by Homan et al. (2018, 20202) and do not need to 
be repeated here. The latter observations (referred to as TMO) 


used an extended antenna configuration and were combined with 
the earlier observations of lower angular resolution by Homan 
et al. (20202) with a low weight of 20% due to fear of possible 
artifacts associated with the lack of short spacings. In the present 
article, in order to better understand their contribution to the 
combined data, we reduced them separately and report about the 
results in the present section. The combined data are obtained 
from Homan et al. (20202) via the CDS (http: / /cdsarc.u-strasbg. 
fr/ viz-bin/cat/J/ A--A/642/493). It should be sufficient to 
summarize their main features in Table 1. 

The 2019 observations were made in two steps, on June 5 
and 14, respectively, with a hybrid antenna configuration 
(C43-10 and C43-9) including ~46 antennas having baseline 
lengths between 83 and 16,000 m. A gap at short spacing 
(~100-200 m), seen in the baseline length distribution of June 
14 (Figure 2(a)), is due to the absence of four antennas (circled 
in red in Figure 2(b)). Four spectral windows were selected to 
observe line emissions of CO(2-1), SiO(5-4), SO2(42 5-3, 3) and 
SCOQ-1). However, the latter two sets of data had insufficient 
sensitivity to contribute significantly new information and their 
analysis is omitted from the present work. 

The data of June 5 and 14 observations of the CO(2-1) and SiO 
(5-4) line emissions have been merged with equal weights and 
grouped in frequency channels corresponding to Doppler velocity 
intervals of 0.95 and 0.34kms~', respectively. Natural weighting 
was used for cleaning, producing maps of 2000 x 2000 pixels 
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Figure 2. (a) distribution of the decimal logarithm of the uv distance (in meters) for the different data sets; (b) antenna configuration (center region); the four antennas 
circled in red were absent for the June 14 observations; (c, d) radial distributions of CO (black) and SiO (red) emissions using TMO (c) and combined data (d). 


covering +4” in both x (pointing east) and y (pointing north), 
with full width at half maximum (FWHM) beam sizes of 31.7 x 
28.0 mas? and 32.4 x 30.0 mas”, respectively. We did not subtract 
the continuum contribution at visibility level, which would have 
been improper when looking over the stellar disk because of 
important absorption. We performed self-calibration but, as it did 
not result in a significant improvement of the quality of the line 
images, the data used here do not include self-calibration. 

The impact of the lack of short spacing below ~100 m was 
found to be different between the CO(2-1) line, with emission 
extending up to some 10" angular distance from the star and 
SiO(5-4) line emission, confined within some 3". Moreover, in 
the case of the CO(2-1) line emission, it was also observed to 
be stronger near systemic Doppler velocity, where the emission 
is dominated by a disk extending to large angular distances. 
The size of the mask, which defines the region in the source 
plane over which the cleaning procedure is applied, must 
therefore be chosen accordingly. We studied in detail images 
obtained with different mask sizes and retained a radius of 2" 
for the SiO(5-4) data and of 1" for the CO(2-1) data, resulting 
in noise levels of 1.34 mJy beam! for SiO(5-4) line emission 
and 1.0mJy beam ! for CO(2-1) line emission. While the 
nominal maximal recoverable scale (MRS) corresponding to 
the C43-10 antenna configuration is 0"46, the impact on 
imaging of the lack of short spacing cannot be simply measured 
by a single number. As the angular distance from the star 
increases, it reduces the detected flux first uniformly, then 
introduces anisotropic distortions of the emission pattern. We 
studied its effect in much detail, both comparing the produced 
images with those given by the combined data and by modeling 
the source as disks of various sizes. While the decline of 
detected flux starts well below the nominal MRS, important 
distortions appear only beyond ~0”6-0”7. For the purpose 
of the present work, we checked carefully that whenever 
TMO data were used, proper account was taken of this issue. 
Channel maps of the SiO(5-4) and CO(2-1) line emissions are 
shown in Figures S1 and S2 of the supplementary material. 


Examples of comparison between TMO and combined data are 
illustrated in Figures 2, 4, 14 and 15. 

We use Cartesian coordinates with the z axis pointing away 
from Earth in addition to the x and y axes in the plane of the sky. 
We refer the Doppler velocity V. to a systemic star velocity of 
—33.6km s! in conformity with earlier work, but its value is 
not known to better than a fraction of a kms ''. We define an 
angular distance from the center of the star as R = J/x? + y? and 
a position angle w measured counter-clockwise from north. 


3. Equatorial Density Enhancement 
3.1. CO(2-1) Emission 


Evidence for the presence of an EDE is known from 
previous work to rest on the observation of the CO(2-1) and 
CO(1-0) line emissions exclusively (Hoai et al. 2019). 
However, a clear distinction between EDE and polar outflows 
could only be made at distances from the star exceeding ~2”. 
The availability of high angular resolution data makes it now 
possible to study the EDE morphology at short distances from 
the star. As illustrated in the upper panels of Figure 3, Doppler 
velocity spectra integrated over successive rings centered on 
the star give clear evidence for the presence of the narrow 
component when R exceeds ~200 mas. It covers Doppler 
velocities between 0 and --2 km s^ !. Remarkably, the neigh- 
boring channels, between —3 and 0 km s | and 2 and 5 km s, 
respectively, are depopulated. This may suggest that very close 
to the star gas is bent toward the equatorial plane, depopulating 
these channels to the benefit of the narrow component. Channel 
maps integrated over the corresponding intervals of Doppler 
velocity are displayed in the lower panels of Figure 3. In the 
blueshifted depopulated Doppler velocity intervals, emission is 
enhanced in the northeastern quadrant. The intervals of high 
|V.| values (—10 to —3 and 5 to 12 kms ) host outflows 
expanding along the polar axes; they are discussed in Section 4. 

In order to better explore the morpho-kinematics of the 
nascent EDE, we display in Figure 4 V. versus w maps of the 
CO(2-1) emission in 100 mas wide rings centered on the star. 
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Figure 3. CO(2-1) line emission, TMO data. Upper panels: Doppler velocity spectra in successive rings, each 100 mas wide, starting from the circle of radius 100 mas 
centered on the star. Lower panels: brightness maps averaged over the velocity intervals associated with the structure of the spectra displayed on the upper panels for 
R > 200 mas, and indicated on top of each panel. 
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Figure 4. CO(2-1) emission: TMO (upper row) and combined (lower row) data. V. vs. w maps in 100 mas wide rings as indicated on top of each panel. The vertical 
lines are at w values of 180°, 270? and 360°. The sine waves shown as reference correspond to pure rotation with a velocity of 4 km s !. Pure expansion of a same 
velocity would be phase-shifted by +90°. The color scales are in units of mJy beam. 


They show that the EDE is already present in the smaller ring, explicitly earlier although Homan et al. (2020a) had remarked 
with R between 100 and 200 mas, and its kinematics is that it is suggested by their SPHERE data. In particular, it raises 
dominated by rotation at a few km s! scale. This evidence for questions in the context of the interpretation given by Homan 
very precocious formation of the EDE had not been made et al. (20202), which would imply an inner EDE radius of ~400 
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Figure 5. CO(2-1) line emission, combined data. Left pair of panels: PV maps in the V. vs. w (panel (a)) and V, vs. R (panel (b)) planes of the projections of the data 
cube defined as 1/2 < R < 5/2 and |V.| < 8 km s. Right pair of panels: distribution in the w vs. R plane of the mean (panel (c)) and rms (panel (d)) values of V. 


over the central Doppler velocity interval, |V.| < 2 km s !. 


mas, corresponding to the orbit of the assumed companion, with 
strong anisotropy. Unfortunately, with presently available data, 
it is difficult to tell apart, in the immediate neighborhood of the 
star, say below R ~ 200 mas, the contributions of expansion, of 
rotation and of phenomena related to the generation of the 
nascent wind: dust formation and shock waves caused by the 
interaction between stellar pulsation and convective cell partition 
(see Section 5). 

Farther away from the star, the EDE is seen to subsist up to 
large distances, as illustrated in the pair of left panels of 
Figure 5 by the PV maps of the line emission seen in the ring 
1"2 « R « 5"2 (see also channel maps of the CO(2-1) 
emission shown in Appendix in both the x versus y plane, 
Figure Al, and the w versus R plane, Figure A2). These maps, 
drawn in the V. versus R and V, versus w planes, respectively, 
show enhanced emission over all values of R and w in the 
central Doppler velocity interval, |V.| < 2 km s !. Similar PV 
maps drawn in small intervals of w and R are shown in 
Appendix (Figures A3 and A4). In each (w, R) bin, with a width 
of 10? x 100 mas, we calculate the mean and rms values of V. 
over the central Doppler velocity interval, |V.| < 2kms !. 
Their values are shown as a function of R and w in the pair of 
right panels of Figure 5. They display remarkable uniformity, 
with the rms value staying near 1 kms! over the whole radial 
range and the mean value near zero. A small w asymmetry of 
the latter, below R ~3”, reveals an inclination of the EDE with 
respect to the plane of the sky consistent with the known 
inclination of the axis of the CSE with respect to the line of 
sight: it reaches values of —.0.3 kms ! in the northwestern 
direction and of ~—0.3kms_! in the southeastern direction. 
This suggests a radial expansion velocity at the scale of 
~0.3/sin10° = 1.7 km s^!, consistent with the low value of 
1.9 km s^! retained in earlier studies (Hoai et al. 2019). The 
rms deviation of 1 kms ! combines the effect of the known 
rms flaring angle of ~15°, sin15? x 1.7 = 0.44 kms | to an 
intrinsic width of ~0.9 km s '. 


3.2. SiO(5-4) Emission 


In order to compare CO(2-1) and SiO(5-4) line emissions in 
the close neighborhood of the star, we display in Figure 6 for 
SiO(5-4) (TMO data) the same distributions as shown in 
Figure 3 for CO(2-1). They show a very different picture than 
those of CO(2-1) emission. The spectra cover a broad range of 
Doppler velocities, independently from the distance to the star, 
with no sign of change beyond 200 mas. The narrow 
component contributes only a very small enhancement. 
Absorption from the outer layer having reached terminal 
velocity at 8-10 km s7" is seen on all four spectra. In contrast 
with CO(2-1) emission, the brightness maps show little 
structure but reveal the presence of a strong depression of 
emission at ~400 mas west from the center of the star. This 
depression has been abundantly studied and commented upon 
in the earlier literature and we report on its properties in 
Section 4, to which we defer possible interpretations. Here we 
simply remark that the different response to temperature of the 
two lines can only have small effects. They differ essentially by 
the value of the energy Ey, of the upper level, 16.60 K for 
CO(2-1) and 31.26 K for SiO(5-4). This quantity enters in the 
expression of the emissivity as a factor e^ ^w/T, where T is the 
temperature. The ratio between CO(2-1) and SiO(5-4) line 
emissions of the values taken by this factor for T — 10, 20, 50 
and >100 K is 4.3, 2.1, 1.3 and «1.2, respectively. 


3.3. Patterns of Enhanced CO(2-1) Emission 


The emission of the CO(2-1) line is known from earlier work 
to display significant fluctuations. Hoai et al. (2019) remarked 
that both CO(1-0) and CO(2-1) emissions in the EDE display 
remarkably similar intensity fluctuations, at the level of 
~+36%, dominated by an eastern enhancement in the form 
of an arc and a southwestern depression, together having the 
appearance of a spiral. Homan et al. (2018) described similar 
data as evidence for a bright nearly face-on spiral structure, 
claiming the first convincing detection of a spiral in the 
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Figure 6. Same as Figure 3 for the SiO(5-4) emission. 


circumstellar environment of an oxygen-rich AGB star, with a 
higher degree of small scale hydrodynamical structure and 
wider spiral arms than in carbon-rich environments. In their 
second paper, Homan et al. (2020a) gave a very detailed 
description of the EDE intensity fluctuations and concluded 
that the observed complex emission pattern strongly resembles 
the typical hydrodynamical gas flow patterns associated with 
the interaction of a mass-losing AGB star with a companion. 
They claimed to have identified the flow along the L1 and L3 
Lagrange points, before they coalesce into the spiral that 
manifests at larger length-scales and argued for the presence of 
a companion located at 0/4 to the west of the star, probably a 
white dwarf, with a mass between 0.65 and 0.8 Mo. 

Channel maps of the combined data of CO(2-1) line 
emission in the twelve frequency channels that bracket the 
systemic velocity of the star (—2.0 < V. < l.8kms ) are 
shown in Figure 7. They display arcs of enhanced emission. 
The complexity of the patterns makes a reliable interpretation 
in terms of well identified individual components difficult. In 
particular it makes it easy to interpret two partly overlapping 
arcs as a single one and an arc that is not centered on the star as 
a spiral arm. The human eye is too subjective a tool for a 
reliable interpretation of the observed patterns. In this context, 
the reality of the spiral discovered by Homan et al. (2018) 
deserves being critically discussed. There is no doubt that 
identifying such a spiral in the pair of central frequency 
channels (centered at 0.05 and —0.27 km gt Figure 7) is a 
natural and reasonable thing to do. However, we note that the 
neighboring channels (centered at 0.37 and —0.59 kms, 


respectively) host patterns suggesting that the source of 
emission is in a plane inclined with respect to the sky plane: 
they shine in the northwestern direction for the redshifted side 
and in the southeastern direction for the blueshifted side, as 
expected for a plane perpendicular to the known axis of the 
CSE. Correcting for this tilt by averaging the brightness over 
the four central channels, (—0.75 « V. « 0.53kms !) we 
obtain the map shown in the left panel of Figure 8. The 
identification of a spiral as proposed by Homan et al. (2018) 
and shown in the central panel of the figure is no longer 
obvious. Assessing the reliability of the identification of such a 
spiral is indeed difficult. As illustrated in the right panel of 
Figure 8, arcs of spirals are excellent approximations to arcs of 
circles having their center offset with respect to the center of 
the star. But the two halves of a single circle correspond to 
different spirals winding in opposite directions; a proper spiral 
winds instead in a single direction and expands indefinitely: an 
arc having a point of closest approach to the star cannot be 
described as a single spiral. 

In summary, the EDE is mostly seen in CO line emission as 
a disk slightly inclined with respect to the plane of the sky, 
formed within some 100 mas from the center of the star and 
having an rms flaring angle of ~15°. Its small inclination 
prevents a precise evaluation of the morpho-kinematics but 
evidence has been obtained for rotation in the central region 
evolving to radial expansion at a low rate of less than 2 km s7! 
beyond a fraction of an arcsec from the center of the star. Such 
a small expansion velocity results probably in part from the low 
penetration of the stellar radiation into the dense disk, limiting 
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clockwise for the blue spiral and counter-clockwise for the red spiral. 


its effectiveness at accelerating dust grains. Assuming a 10” 
disk radius (Figure 6 of Homan et al. 2018), it implies that the 
disk is at least thirty centuries old. Arcs of enhanced emission 
are observed on the disk but their interpretation as forming a 
spiral has been disputed. 


4. Polar Outflows 
4.1. CO(2-1) Emission 


The polar outflows cover a broad region of the sky up to 
large velocities. This is illustrated in Figure 9, which displays 
V. versus R maps integrated over position angles for the 
blueshifted and redshifted hemispheres separately. A simplified 


schematic of the morpho-kinematics of the CSE is shown in the 
right panel of Figure 9. 

The polar outflows were observed earlier to display density 
fluctuations at the level of ~+26%, and to reveal both an 
asymmetry between blue- and redshifted hemispheres and 
significant differences between CO(1-0) and CO(2-1) emis- 
sions (Hoai et al. 2019). The higher angular resolution of the 
present data allows for an improved study of the morphology of 
the enhancements of CO(2-1) line emission. We studied in 
great detail the data cube (combined data) covering 51 
frequency channels (Doppler velocities between —8.2 and 
--8.2kms !) and projected angular distances from the star 
between 172 and 5"2 (Figures Al to A4 of the Appendix). In 
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some cases where one might have feared that the steep radial 
decline of the intensity would have been a source of bias, we 
examined the brightness divided by a power fit to its radial 
distribution in the relevant frequency channel, found to have an 
index between —0.2 and —0.4. In no case was it considered to 
be a source of significant bias. 

Ideally, one would like to describe each channel map as a set 
of well-identified separated arcs of enhanced emission. We 
attempted to do so by using an algorithm lumping together the 
brighter pixels. However, we failed to obtain a satisfactory result, 
the complexity of the patterns preventing reliable pattern 
recognition. Instead, for each frequency channel, we inspected 
carefully the different projections of the data cube and tried to 
describe the observed patterns as simply and reliably as possible. 
Figure 10 displays a pair of representative examples; they have 
not been selected to show particularly clear patterns. There are 
many examples of clearer cases as well as of more obscure cases. 

Rather than giving a detailed review of our observations, we 
simply list the main statements which we can make with 
confidence: in each map, the pattern of enhanced emission can 
be described as a set of arcs; they are too numerous, and 
frequently overlapping, to allow for reliable pattern recogni- 
tion; they can be approximated by arcs of circles covering a 
few tens of degrees, typically less than 180°; the contrast is 
small, with a peak-to-valley ratio not exceeding a factor of 2; 
the arcs always wind around the star, implying that the center 
of the associated circle is offset by less than a radius, typically 
half a radius, from the center of the star; their radii cover a 
broad range of values, from as small as permitted by the 
angular resolution to as large as permitted by the size of the 
maps under study, namely from a fraction of an arcsecond up to 
over 5”; very similar patterns are observed in adjacent 
frequency channels, covering ranges up to a few km s !. These 
comments complement usefully those made in this context by 
Homan et al. (20202). 


Finally, we remark that the V. versus R map shown in 
Figure 5(b) (see also Figure A3 of the Appendix) provides 
precious information on the dependence of the radial expansion 
velocity on stellar latitude a. It was assumed to be of the form 
V = Veg + (Vole — V) X sin? a in Hoai et al. (2019), where 
Veg and Vpole are the expansion velocities at equator and pole, 
respectively. As illustrated in Figure 11, a small volume of gas 
at distance r from the star, at latitude a and having a radial 
velocity V has coordinates V, — Vsina and R=rcosa. 
Namely, a spherical shell expanding at radial velocity Vo and 
observed when it reaches a radius Ro is seen as a circle of unit 
radius in the PV map V,/Vo versus R/Ro. In contrast, a ring 
expanding at the equator at velocity V; is seen as a point on the 
R axis at V; / Vo. Understanding the transition between the polar 
outflows and the EDE implies understanding the transition 
between these two regimes. Namely, a shell emitted at some 
time with velocity Vo at the pole and velocity V; at the equator 
is seen as a curve of equation V/V, = Vi/Vo + (1 — Vi/Vo)yf (8) 
where J is the polar angle, with f(0)=0 and f(90°)=1. 
However, ) is equal to the stellar latitude o because the 
reference radius Ro used to define the axis of abscissa is the 
distance covered at velocity Vo, used as reference velocity on 
the axis of ordinate. Namely the curve representing the shell in 
the PV diagram is a direct representation of the function f(a) 
that describes the dependence on stellar latitude of the 
expansion velocity. The left panel of Figure 11 shows three 
examples of function f(o), all three of the form sin"o, with 
n= > 1 and 2 respectively (in addition to the trivial cases 
f(a) =0 and 1). Comparison with low latitude observations, as 
shown in the right panels of Figure 11 (see also Figure S3), 
favors clearly values of n smaller than 1, in contradiction with 
the arbitrary assumption made in Hoai et al. (2019). 
Specifically, it suggests that the outflows are close to portions 
of spherical shells when emitted at large stellar latitudes, and 
that their radial velocity decreases rapidly at small latitudes 
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Figure 11. Latitudinal dependence of the radial expansion velocity. Left: schematics. Center and right: two examples of V. vs. R maps (CO(2-1) combined data) at 
w = 65° and 205°, respectively (from Figure S3). The red lines in the central and right panels guide the eye on the flows of maximal emission. 


where their interaction with the disk becomes important. This is 
important information for modeling the hydrodynamical flow, 
however a task beyond the scope of the present work. 
Tentatively, the properties of the enhanced patterns of CO 
line emission in the polar outflows suggest that they take the 
form of shells of episodic mass ejections and associated shock 
waves triggered by pulsations of amplitudes enhanced on the 
surface of convective cells. The confinement of such cells 
within solid angles at steradian scale limits the patterns 
observed in the channel maps to finite arcs rather than full 
circles. Their radial velocity is maximal at the center of the cell 
and drops near its edges, causing the intersection of the shell 
with the sky plane of a given frequency channel to appear 
eccentric. Moreover, their mean radial velocity increases from 
the equator to the poles, implying a similar decrease of their 
density and causing the polar depressions described by Hoai 
et al. (2019) and the appellation of “bi-conical outflows” 
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adopted by Homan et al. (2020a). Their patterns are best 
revealed when looking at the V. versus R distribution of the line 
emission in a wedge of position angles not exceeding the size 
of the convective cell (as done in the pair of right panels of 
Figure 11); integrating over the whole range of position angles 
or over a large radial distance smoothens out the fluctuations 
resulting from the episodic nature of the mass ejections, which, 
as we shall see later, succeed each other at time intervals at the 
scale of a few decades. At the present stage, such an 
interpretation is highly speculative but we shall see that further 
observations described in the reminder of the present section 
will confirm its likely validity. 


4.2. SiO(5-4) Emission: General Features 


In contrast to CO(2-1) emission, S1O(5-4) emission is confined 
within radial distances of ~3”, as shown in Figure 12, which 
displays radial distributions for both SiO(5-4) and CO(2-1) 
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Figure 13. SiO(5-4) emission (combined data) in the redshifted (upper row) and blueshifted (lower row) hemispheres and the radial interval 05 < R < 23. From left 
to right: Doppler velocity spectrum, S(V-); radial distribution, U(R); ratio of the observed intensity to U(R) x S(V-.), properly normalized, in the x vs. y, V. vs. R and V, 


vs. w planes, respectively. 


emissions. The emissions of the two lines are indeed expected to 
differ in several respects: the SiO emission decreases with 
distance to the star, both because of condensation on dust grains 
and of dissociation by UV interstellar radiation; moreover, it is 
an efficient shock tracer. 

As already remarked, and clearly apparent on the distribu- 
tions of SiO(5-4) emission displayed in Figure 13, one observes 
a significant asymmetry between the blueshifted and redshifted 
hemispheres. In particular the former reveals the presence of a 
strong depression of emission at some 400—500 mas west of the 
center of the star that is challenging interpretation. In a first 
step, we limit therefore our considerations to the morpho- 
kinematics of the red hemisphere. To illustrate the small 
inhomogeneity of the data cube, we show in Figure 13 the 
radial distribution U(R), averaged over position angle and 
integrated over V., and the Doppler velocity spectrum S(V.), of 
the volume of the data cube defined as 0"5 < R < 2”3,0 < V. 
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< 8kms !; we then compare the observed intensity to the 
product U(R) x S(V.), properly normalized, in the right panels 
of the same figure. At the same time as the result reveals clearly 
the morpho-kinematics of the observed deviations with respect 
to a smooth distribution, it shows that their amplitude is 
relatively modest, typically at the level of +20%. Both U(R) 
and S(V.) show no significant fluctuation and the lumps of 
enhanced emissions have typical sizes at the scale of a fraction 
of an arcsec. In contrast, on the V. versus R and V, versus w 
maps, they cover several km s !. 

Repeating the same analysis for the blueshifted hemisphere, 
we obtain the result illustrated in the lower row of Figure 13. 
The main difference with the redshifted hemisphere is a lower 
intensity, by ~33%, dominated by an important depression in 
the western quadrant causing a deviation twice as large with 
respect to the smooth distribution. If it were not for this 
depression, in the approximation of neglecting fluctuations at 


Research in Astronomy and Astrophysics, 24:035009 (23pp), 2024 March 


R-(0.1,0.2) 


R-(0.2,0.3) 


200 400 
to (degree) 


R=(0.3,0.4) 


4 
2 
0 


Nhung et al. 


R=(0.4,0.5) 


R=(0.5,0.6) 


AO =- NW A 


200 400 


Figure 14. SiO(5-4) TMO (upper row) and combined (lower row) data. V, vs. w maps in 100 mas wide rings as indicated on top of each panel. Note that the w axis is 
offset by 90°. A triangle and an ellipse, corresponding approximately to the location of the voids of emission, are shown as reference in all panels. The color scale is in 


mJy beam™!. 


the level of +20%, both hemispheres would reveal similar 
morpho-kinematics, which could be described by a simple 
nearly spherical model as previously remarked by Homan et al. 
(2020a). It is therefore important to study it in detail, which we 
do in the following sub-section. 

Channel maps of the SiO(5-4) line emission are shown as 
supplementary material (Figure S3). Also shown are PV maps 
in the V, versus R plane (Figure S4) and in the V, versus w 
plane (Figure S5), which display episodic and lumpy emission. 


4.3. SiO(5-4) Emission in the Close Neighborhood of 
the Star 


Qualitatively, the upper and lower rows of Figure 13 display 
some similarities: below R ~1”, they both show enhanced 
emission in the polar outflows, at V, ~7 and —5 km s 
respectively, with the former covering the southern quadrant 
and the latter covering a broader interval of position angles, 
somewhat enhanced in the northeastern hemisphere. In 
contrast, beyond R ~1”, emission is enhanced near the EDE. 

Figure 14 displays V. versus w maps of SiO(5-4) emission in 
100 mas wide rings centered on the star, in the same format as 
was shown for CO(2-1) emission in Figure 4. Comparing the 
TMO and combined data images, and accounting for the 
different angular resolutions, they clearly reveal the presence in 
the redshifted hemisphere of an outflow covering the southern 
quadrant and of a depression covering the northern quadrant. 
The picture is less clear in the blueshifted hemisphere where 
enhanced emission is seen to cover velocity and position angle 
values closer to those displayed by the disk, while the western 
void shows up clearly in both TMO and combined data images. 

While qualitative similarities between both outflows and 
associated voids are remarkable, as clearly illustrated in 
Figure 15, important differences must be noted: (1) the western 
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blueshifted depression of S1O(5-4) emission is stronger than the 
northern redshifted one, with a peak-to-valley ratio of ~6 
compared with ~3; (ii) the northern SiO(5-4) depression starts 
earlier than the western one and is already present in the 
interval 071 < R < 0"2 where the western depression is not 
only absent but even partly replaced by a clear enhancement; 
ii) the southern outflow, covering large Doppler velocities, is 
emitted in the redshifted hemisphere close to the line of sight 
while the blueshifted outflow, covering small Doppler 
velocities, is likely to interact with the nascent EDE in the 
equatorial plane. Indeed, Figure 16 shows that the CO(2-1) 
emission of the EDE displays a western void of emission 
clearly associated with the edge of the blueshifted void of SiO 
(5-4) emission, as had been previously remarked by Homan 
et al. (20202). 

These observations underscore the complexity of the 
dynamics at stake and suggest that the interaction between the 
blueshifted outflow and the nascent EDE may play an important 
role. Particularly puzzling are the properties displayed by the 
blueshifted western void of SiO(5-4) emission. 

In order to explore these in more detail, we take advantage of 
its confinement to a relatively narrow interval of position angle, 
240? < w < 300° to study the distributions of the CO and SiO 
line emissions in the V. versus R plane, integrated over w, 
which are displayed in the left pair of panels of Figure 17. Both 
give evidence for a mass ejection in the form of a shell having a 
mean radius between 0"5 and 1" and a mean radial velocity 
between 5 and 10 kms |. The void of emission, for both SiO 
and CO lines, is clearly seen in its wake, suggesting that it 
simply results from the radial velocity of the shell significantly 
exceeding that of the wind that followed it. As a distance of 1" 
is covered in the plane of the sky in ~60 yr at a velocity of 
10kms_', the mass ejection occurred a few decades ago. We 
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Figure 16. EDE region, combined data. CO(2-1) (upper row) and SiO(5-4) (lower row) emissions: channel maps centered at the Doppler velocities indicated on top of 
each panel. 
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Figure 17. Distributions of the line emission of SiO(5-4) and CO(2-1) over the V. vs. R plane, integrated over position angles in the interval 240° < w < 300° (left pair 
of panels) and 280? < w « 360? (right pair of panels). In each pair of panels SiO is left and CO is right. 


note that the distribution of the CO line emission suggests the SiO gas has significantly decreased at the corresponding radial 
presence of an earlier mass ejection that had taken place again a distance. In line with the interpretation given at the end of 
few decades before the more recent one. It is not expected to be Section 4.1 of the patterns observed in the CO line emission of 
visible on the distribution of SiO line emission as the density of the polar outflows, Figure 17 suggests the occurrence of mass 
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ejections triggered by pulsations of amplitudes enhanced both 
episodically, with a timescale of a few decades, and in solid 
angle, at the scale of a steradian corresponding to the presence 
of a large convective cell. The angular resolution of the present 
observations cannot resolve the fine time structure associated 
with the pulsation period, ~110 days. 

A significant difference between the distributions of SiO and 
CO line emissions displayed in Figure 17 is the evidence for 
SiO emission to have resumed after the mass ejection, covering 
some 100-200 mas from the center of the star, while CO 
emission has not. Indeed, we saw in Figure 3 that close to the 
star and in the Doppler velocity interval —3 <V, «0 kms! 
CO line emission is instead enhanced in the northeastern 
quadrant. In contrast, we saw from Figure 6 that in the same 
region of the data cube, SiO emission covers the whole range 
of position angles. Understanding the reason for such a 
difference would require understanding the mechanism that 
governs the formation of the equatorial disk as well as 
observations of higher sensitivity and of other line emissions, 
both of which are lacking. We note that the ejected shell covers 
latitudes reaching down to the disk, implying that it interacts 
with it, making it even more difficult to provide a reliable 
interpretation of the mechanism at stake. 

With the aim of exploring the properties of the northern 
depression, we repeated the analysis made for the western void 
in the interval of position angles 280? <w « 360°. The result is 
illustrated in the right pair of panels of Figure 17. Here again, 
in the range of large Doppler velocities, V. 7 kms ', the 
depression is seen in the wake of a shell-like enhancement of 
emission, much clearer on the CO distribution than it is on the 
SiO distribution, and an earlier shell is clearly visible on the 
CO distribution. Also, as for the western void, SiO line 
emission is seen to have resumed close to the star after the mass 
ejection, but CO line emission has not. The features displayed 
by the northern depression are qualitatively similar, but much 
less clear, than they are by the western void. 

While appealing, the interpretation of the western void 
proposed in the present section leaves a number of questions 
unanswered, in particular concerning the mechanism respon- 
sible for the formation of the equatorial disk. The complexity of 
the observed features prevents claiming its validity and unicity 
with strong confidence. However, it questions the validity of 
the interpretation proposed by Homan et al. (20202) in terms of 
photodissociation by a white dwarf companion: such a 
companion should have V. close to zero, like the disk that it 
is meant to have produced, and should irradiate both the 
blueshifted and redshifted hemispheres. Additional observa- 
tions of high sensitivity and angular resolution, also including 
other line emissions, are required in order to ascertain the 
validity of either of the two scenarios. 

Additional figures of relevance to the description of the 
region of the western depression are displayed in Figures S7 to 
S9 of the supplementary material. 
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4.4. SiO(5-4) Emission: Patterns of Enhanced Emission 


In order to compare CO(2-1) and SiO(5-4) line emissions, 
we must limit the radial range and compensate for the different 
radial declines. In practice, we consider the ring 
0"7 < R < 2”0, also excluding the western sextant 240? < 
w « 300° in the blueshifted hemisphere, and divide the 
observed intensity by a fit of the R distribution to a form AR" as 
was done in Section 4.1. Examples of the resulting maps are 
shown in Figure 18, one in the blueshifted hemisphere, with 
b — —0.20 for CO(2-1) and —1.52 for S1O(5-4), and the other 
in the redshifted hemisphere, with b — —0.76 for CO(2-1) and 
—1.20 for SiO(5-4). The normalized brightness has a mean of 
unity by construction and an rms width that provides a measure 
of the intensity contrast, larger for CO(2-1) than for SiO(5-4): 
0.30 and 0.31 for the former, 0.15 and 0.22 for the latter, in the 
redshifted and blueshifted hemispheres, respectively. The 
correlation between the CO(2-1) and SiO(5-4) normalized 
maps is illustrated in Figure 18 by plotting, for each pixel, the 
intensity of the latter versus that of the former. A fit of the form 
Fon (S10) = Fo4-Go x Fio (CO) gives a normalized x. with 
the values of Je in the examples shown in Figure 18 being 
0.14 and 0.21, in the red- and blueshifted hemispheres, 
respectively; the value of Go, which measures the strength of 
the correlation, is 0.18 in both hemispheres. The pair of 
examples displayed in Figure 18 calls for a number of remarks 
of general validity: the maps of CO(2-1) and SiO(5-4) line 
emissions display more similarities in the redshifted hemi- 
sphere than in the blueshifted hemisphere; small enhancements 
of emission are seen on both maps, less contrasted on the SiO 
map than on the CO map; some are common to both maps, 
others are seen on only one of the two maps; in contrast to the 
higher |V.| range illustrated in Figure 10, clear arcs are rarely 


z 


observed. 


5. High Doppler Velocities 


For a detailed study of the neighborhood of the star and of 
the nature of the high velocity wings, we use the high 
resolution TMO data in the form of Doppler velocity spectra of 
SiO(5-4) emission, integrated in 20 mas broad annular rings 
centered on the star with mean radii increasing from 10 to 
390 mas in steps of 20 mas. On each of these we define values 
of the full width measured at 10% of maximum, FWTM, as 
sketched in Figure 19(a). This choice of 10% is a compromise 
between minimizing systematic uncertainties and staying well 
above noise and its precise value is of no relevance to the 
validity of the result. The dependence of FWTM on R gives 
evidence for confinement within some 200 mas (Figure 19(b)). 
It is well described by a fit of the form 21.3+19.0 exp(— 
R/0.10) km s '. The first term corresponds to the velocity 
range covered by the bulk of the CSE, the second term to a 
region of high velocities just above the photosphere, as 
observed in other oxygen-rich AGB stars, and commonly 
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Figure 18. Comparing CO(2-1) and SiO(5-4) emissions (combined data) in the polar outflows: two examples (in the redshifted hemisphere, upper row and in the 
blueshifted hemisphere, lower row). Left pair of panels: normalized channel maps for V. = 3.6 km s ! and V, = —3.5 km s !. We have circled in black some features 
common to both the CO and SiO maps and in red others that are seen in only one of the two maps. Central panels: normalized R and w distributions for CO(2-1) 
(black) and SiO(5-4) (red), respectively. Right: correlation between the normalized brightness of the two lines; the line is the best fit result (see text). 
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Figure 19. (a): definition of the quantity FWTM (see text) taking as example the Doppler velocity spectrum of SiO(5-4) line emission measured in the ring 
60 < R < 80 mas (see Figure S6 of the supplementary material). The blue line is used to extrapolate the spectrum (using a linear ordinate scale), its intersection with 
continuum level defining Vmax (see text). (b): dependence of FWTM on the ring radius, R. The line is a fit 21.3+19.0 exp(— R/0.10) km s. (c): comparison between 
the measured dependence on R of the mean value of Vnax and the best fit F(R) = (1.13—0.11R). (d): comparison between the measured dependence on w of the mean 


value of Vmax and the best fit Vo(w) = [11.14-0.8 sin(w —71?)]. 


interpreted as resulting from shock waves produced by the 
interaction of stellar pulsation and convective cell partition (see 
for example Hoai et al. 2020 and Nhung et al. 2022). We also 
estimate the maximal Doppler velocity Vmax, reached in the 
redshifted hemisphere by the bulk of the expanding CSE, by 
extrapolating the upper end of the spectrum down to continuum 
level (essentially zero) using the combined data. The depend- 
ence of Vmax on R and w is well described by a fit of the form 
Vinax = (1.13—0.11R) x [11.14-0.8 sin(w—71°)]. The rms devia- 
tion between data and fit is 0.36 kms '. Panels c and d of 
Figure 19 compare these two factors to the measured mean 
values of Vinax, averaged over w and R, respectively. The 
dependence on w corresponds well to what is expected from the 
tilt of the CSE axis, pointing ~20° east of south in the red 
hemisphere. 

This analysis gives therefore evidence for two major 
contributions to the high velocity wings of the SiO(5-4) Doppler 
velocity spectrum: one, extending up to some 12kms |, is 
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simply the upper edge of the velocity distribution of the polar 
outflows; the other, reaching up to 20 km s^, is confined to the 
close neighborhood of the star and corresponds to what is 
observed in most other oxygen-rich AGB stars where this region 
of the CSE has been carefully explored. The interpretation given 
by Tuan-Anh et al. (2019) of these high velocity wings, using 
data of lower angular resolution, was mixing both contributions. 
This is clearly illustrated in Figure 20, which displays PV maps 
of |V.| versus R for both line emissions and both hemispheres 
separately. In all four cases there is evidence for a clear 
separation between the two contributions, even if significant 
differences exist between the CO(2-1) and SiO(5-4) lines as well 
as between the redshifted and blueshifted hemispheres. Another 
illustration is given by the intensity maps displayed in the left 
pair of panels of Figure 21 for 8 < |V.| < 12 kms | (combined 
data) and in the right pair of panels of the same figure for 12 < 
|V.| < 20kms~' (TMO data). The former gives evidence for the 
above-mentioned asymmetry between blueshifted and redshifted 
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Figure 20. High Doppler velocities (combined data): PV maps of |V.| vs. R for both hemispheres and both lines separately, as indicated on top of each panel. 
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Figure 21. Left: intensity maps of the SiO(5-4) emission (combined data) measured in the Doppler velocity intervals 8 < |V.| < 12 km s! 
the SiO(5-4) emission (TMO data) measured in the Doppler velocity intervals 12 < |V.| < 20 km s7 


hemispheres along the projection of the CSE axis, while the 
latter displays instead a small elongation in the perpendicular 
direction where they are better distinguishable from the bulk. 
This latter high V. component seems however different from 
what we know of other stars by the absence of redshifted 
absorption, usually interpreted as evidence for in-falling gas, 
over the stellar disk (Figure 22). Unfortunately, the lack of 
reliable light curve measurements prevents commenting on a 
possible relation with the stellar phase at the time of 
observations. 


6. Summary and Discussion 


Collecting the information obtained in both the present work 
and earlier publications to draw a credible picture of the 
morpho-kinematics of the CSE of EP Aqr is challenging. The 
interpretation given by Homan et al. (20202), giving a major 
role to a white dwarf companion currently located at ~400 mas 
west of the star, has been shown to be not unique and, 
accordingly, to be less reliable than previously thought, by the 
results of the analyses made in the present work: the evidence 
for a spiral enhancement of CO(2-1) emission and the 
interpretation of the lumps of CO(2-1) emission in the 
neighborhood of the star illustrated in Figure A12 of Homan 
et al. (20202) have been shown to become invalid when 
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considering a slightly broader Doppler velocity interval. The 
detailed study of the blueshifted western depression of SiO(5- 
4) emission presented in Section 4.3 failed to support an 
interpretation in terms of dissociation by UV radiation from a 
white dwarf companion. This does not mean that one can 
exclude the presence of a companion and its possible role in 
shaping the CSE. In particular, there is no reason to doubt that 
EP Aqr, as most other stars, has several planets orbiting it, 
some of which may already have been engulfed. However, the 
complexity of the observed morpho-kinematics, and the 
resulting difficulty to give a reliable interpretation of its many 
features, prevent claiming a credible model of the mechanisms 
at stake with sufficient confidence. 

In general, the global morpho-kinematics of CO(2-1) 
emission observed in the present work confirms the results of 
earlier analyses. The evidence for approximate axi-symmetry 
about an axis making an angle of ~10° with the line of sight 
and projecting ~20° west of north on the plane of the sky has 
been strengthened by observations of both the EDE and polar 
outflows. The latitudinal dependence of the expansion velocity, 
evolving from ~2 to ~10kms~! from equator to poles, has 
been shown to imply a sin"o term with n smaller than unity, in 
contrast with the assumption of n —2 made by Hoai et al. 
(2019) (Section 4.1). Evidence for enhancements of emission, 
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Figure 22. Doppler velocity spectra of SiO(5-4) (upper panels) and CO(2-1) (lower panels) line emissions (TMO data) integrated in rings as indicated on top of the 


upper panels. 


with a typical peak-to-valley ratio below a factor 2, has been 
found in both EDE and polar outflows, as well as in both CO(2-1) 
and SiO(5-4) line emissions, often, but not always, in 
coincidence. Their morphology in the CO line emission of the 
polar outflows, taking the form of eccentric arcs, suggests that 
they are caused by shells of mass ejections covering solid angles 
at steradian scale. 

The birth of the equatorial enhancement observed in CO(2-1) 
emission has been shown to occur very close to the star, first 
dominated by rotation, then by expansion, becoming particu- 
larly clear at distances exceeding ~200 mas below which 
important effective line broadening and dust formation are 
known to be present (Section 3.1). Lack of evidence for a 
similar enhancement in SiO(5-4) emission is puzzling; only a 
very small enhancement is visible at distances exceeding 
7-300 mas (Section 3.2). This suggests that a disk is formed at 
short distances from the star probably caused by the presence or 
the recent engulfment of a companion (Privitera et al. 2016); 
but this is pure speculation. An intrinsic rms width of nearly 
1 kms! has been found to contribute to the EDE flaring of 
~40° FWHM evaluated earlier (Hoai et al. 2019). 

The strong difference between the morpho-kinematics displayed 
by the SiO(5-4) and CO(Q-1) emissions (Section 4.2) was 
understood as being largely due to the different radial ranges 
being explored, the former declining rapidly as a result of accretion 
by dust grains and dissociation by the interstellar UV radiation. In 
a first approximation, SiO(5-4) emission can qualitatively be 
described by the standard mass-loss mechanism at play in other 
oxygen-rich AGB stars (Hófner & Freytag 2019), with a typical 
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terminal radial expansion velocity at the scale of l0kms ! and 
some prolateness caused by the presence of the EDE. 

Indeed, a careful inspection of the morpho-kinematics of the 
S1O(5-4) line emission in the close neighborhood of the star has 
revealed the presence of two distinct broad outflows hosting 
episodic mass ejections at a timescale of a few decades, 
covering solid angles at steradian scale, probably associated 
with shock waves and leaving a depression in their wake: if it 
were not for the presence of the EDE, the morpho-kinematics 
of the CSE could be understood as resulting from the standard 
mechanism of wind generation described by current state-of- 
the-art hydrodynamical models (Freytag & Hófner 2023). 
However, the presence of the EDE causes the CSE to take the 
form of two polar outflows having an increasing radial velocity, 
and an accordingly decreasing density, when spanning from the 
equator to the poles. 

Such a model, suggested by the results of the present 
analyses, is at significant variance with that proposed by 
Homan et al. (20203). However, in order to claim its 
superiority, we would need to understand the details of what 
precisely causes the progressive disappearance of SiO 
molecules, of what governs the morpho-kinematics of the 
EDE, of what kind of dust takes part in the acceleration 
process, of what are the relative roles of stellar pulsation and 
convective cell partition in giving the initial boost to the wind, 
and of which kind of pattern of shock waves is present. 

Only then could one decide with confidence which of the two 
models is closer to reality. In the present stage, we need 
additional ALMA observations of high angular resolution and 


Research in Astronomy and Astrophysics, 24:035009 (23pp), 2024 March 


sensitivity to tell them apart. Particularly informative would be 
observations of the emission of molecular lines that probe 
preferentially the very close environment of the star, such as 
from vibrational excited states. New infrared observations of the 
dust, together with simultaneous monitoring of the light curve, 
are required to possibly reveal correlation with the stellar phase. 
Finally, the presence of high velocity wings, extending up to 
~20kms', previously noted in earlier work, has been 
confirmed and shown to include two distinct components 
(Section 5): upper velocity ends of the polar winds, reaching 
typically 12 km s^! but with weak tails reaching well above; and 
effective line broadening confined to distances shorter than 
200 mas from the center of the star. The earlier interpretation, 
given by Tuan-Anh et al. (2019), was failing to make the 
distinction between these two components. The observed line 
broadening, consistent with similar observations in all other low 
mass-loss rate oxygen-rich AGB stars observed with sufficient 
angular resolution and sensitivity, is believed to be due to lumpy 
emission of shock waves associated with convective cell 
partition and stellar pulsation, but a detailed picture is lacking. 
In general, presently available observations do not allow for 
telling apart, in the close neighborhood of the star, the 
contributions of pulsations, convective cell partition, EDE 
formation and dust condensation. EP Aqr is clearly a star from 
which one should learn more than currently feasible; it deserves 
being studied in as many details as possible with state-of-the-art 
instruments; it still leaves us with many open questions, the 
answers to which would be of great help to improve our 
understanding of the mechanisms at stake in the generation and 
development of the wind of oxygen-rich AGB stars. 
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Appendix 


We show below channel maps of the CO(2-1) emission in 
both the x versus y plane (Figure A1) and the w versus R plane 
(Figure A2). Also shown are PV maps in the V. versus R (Figure 
A3) and the V. versus w planes (Figure A4). 
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Figure A1. Channel maps of CO(2-1) emission (combined data). Central Doppler velocities (km s~') are indicated for each channel on the upper-left corner of the 
panel. 
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Figure A2. Channel maps of CO(2-1) emission (combined data) in the w vs. R plane. Central Doppler velocities (km s~') are indicated for each channel on top of the 
panel. 
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Figure A3. Projections of the data cube of CO(2-1) emission (combined data) in 36 intervals of w (the central value in degree is indicated on top of each panel) on the 
V. vs. R plane (V.| < 8kms !, 1.2 < R < 5/2). 
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Figure A4. Projections of the data cube of CO(2-1) emission (combined data) in 40 intervals of R (the central value in arcsec is indicated on top of each panel) on the 
V. vs. w plane (|V.| < 8 km s !, 0 < w < 360°). 
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